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Abstract 

The quasinormal modes (QNMs) of Dirac field perturbations of a Reissner-Nordstrom black hole in 
an asymptotically Anti-de Sitter spacetime are investigated. We find that both the real and imaginary 
parts of the fundamental quasinormal frequencies for large black holes are the linear functions of 
the Hawking temperature, and the slope of the lines for the real parts decreases while that for the 
magnitude of the imaginary parts increases as the black hole charge increases. According to the 
Anti-de Sitter /Conformal Field Theory (AdS/CFT) correspondence, the fact shows that different 
charge presents different time scale in three-dimensional CFT. Another interesting result is that 
the quasinormal frequencies become evenly spaced for high overtone number, and in the spacing 
expressions the real part decreases while the magnitude of the imaginary part increases as the charge 
increases. We also study the relation between quasinormal frequencies and angular quantum number 
and find that the real part increases while the magnitude of the imaginary part decreases as the 
angular quantum number increases. 
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I. INTRODUCTION 



AdS black holes have a direct inter- 
According to the AdS/CFT 



It is well known that QNMs possess discrete spectra of complex characteristic frequencies 
which are entirely fixed by the structure of the background spacetime and are irrelevant of 
the initial perturbations 0. Thus, it is believed that one can directly identify a black hole 
existence by comparing QNMs with the gravitational waves observed in the universe, as well as 
test the stability of the event horizon against small perturbations. Meanwhile, it is generally 
believed that the study of QNMs may lead to a deeper understanding of the thermodynamic 
properties of black holes in loop quantum gravity since the real part of quasinormal 
frequencies with a large imaginary part for the scalar field in the Schwarzschild black hole is 
equal to the Barbero-Immirzi parameter 00 00, a factor introduced by hand in order that 
the loop quantum gravity reproduces correct entropy of the black hole. 

Motivated by the AdS/CFT correspondence, the study of QNMs in AdS spacetime becomes 
appealing recently. It was argued that the QNMs o: 
pretation in terms of the dual conformal field theory 
correspondence, a large static black hole in asymptotically AdS spacetime corresponds to an 
(approximately) thermal state in CFT, and the decay of the test field in the black hole space- 
time corresponds to the decay of the perturbed state in CFT. The dynamical timescale for the 
return to the thermal equilibrium can be done in AdS spacetime, and then translated onto the 
CFT, using the AdS/CFT correspondence. The QNMs in AdS spacetime were first computed 
for a conformally invariant scalar field by Chan and Mann [lflj . Subsequently, Horowitz and 
Hubeny 0] presented a numerical method to evaluate the quasinormal frequencies directly and 
made a systematic investigation of QNMs for the scalar field in the Schwarzschild AdS black 
hole. Then, many authors studied QNMs for the scalar, electromagnetic and gravitational 
perturbations in various asymptotically AdS black holes -24|. Although QNMs for the 
scalar, electromagnetic and gravitational field perturbations in the Schwarzschild AdS and the 
Reissner-Nordstrom AdS black-hole backgrounds have been investigated extensively, QNMs for 
;he Dirac field perturbations in the Schwarzschild AdS black hole were first studied in Refs. 
25I ] recently. Considering the Reissner-Nordstrom AdS (RNAdS) black hole spacetime pro- 
vides a better background than the Schwarzschild AdS geometry, in this paper we will extend 
the study in Ref. |2y] to the RNAdS black hole and see how the black hole charge affects the 
Dirac QNMs. 

The organization of this paper is as follows. In Sec. 2 the decoupled Dirac equations and the 
corresponding wave equations in the RNAdS spacetime are obtained by using Newman-Penrose 
formalism. In Sec. 3 the numerical approach to compute the Dirac QNMs is introduced. In 
Sec. 4 the numerical results for the Dirac QNMs in the RNAdS black hole are presented. The 
last section is devoted to a summary. 
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II. DIRAC EQUATIONS IN THE REISSNER-NORDSTROM ANTI-DE SITTER 
SPACETIME 



The line element of the RNAdS black hole can be expressed as 



with 



ds 2 = fdt 2 - jdr 2 - r 2 (d6 2 + sin 2 6d V 2 ), (2.1) 



/ = 1 — + ~T + ~E2> (2.2) 
r r z R l 



where M, Q and R represent the mass, charge and anti-de Sitter radius, respectively. Hereafter 
we will take R = 1 and measure everything in terms of R. The spacetime causal structure 
depends on the zeros of /. Changing the parameters M and Q, the function / may have none, 
one or two positive zeros. In general case, / has two simple real, positive roots r_ and r + , 
but in the so-called extreme case, / has one double positive zero r + (because r_ = r+). The 
horizons r_ and r + with r_ < r + , are called Cauchy and event horizons respectively. The 
Hawking temperature of the black hole is given by 

T H = ± = !~ W + fri , (2 .3) 
2tt 47rr + 

where k is the surface gravity. The mass parameter M is related to the charge Q and the event 
horizon radius r + by the relation 

1 (. , .3 , Q 2 



M=-i [ r + + ri + f-j } (2.4) 

and the extremal value of the black hole charge, Q ex t, is given by 

Qlxt = r 2 + (l + 3r 2 + ) ■ (2-5) 
In a curved spacetime the Dirac equations j22j are described by 

V2V BB ,P B + ifiQ B , = 0, 

y/2V bb>Q B + wPb> = 0, (2.6) 

where V bb< is covariant differentiation, P B and Q B are the two-component spinors representing 
the wave function, Pb> is the complex conjugate of P B , and \i is the particle mass. In order 



to separate t 
formulism 128 



re Dirac equations in the RNAdS spacetime (|2.1|) by using Newman-Penrose 
we take the null tetrad as 

2 

P = (^, 1, 0, 0), 



1/ A 

2< X ' —r 



i, o, 0) 



0,0,1,-^), (2.7) 
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with 



A = r 2 - 2Mr + Q 2 + r 4 



(2.8) 



and set the two-component spinors as 



P° = -R_ 1/2 SL 1/2 e 
P 1 - 

Q 1 ' 

~0' 



i(u}t—mip) 

/ 2 S +1/2 e-^-"^, 



Q 



1 



-i(ujt— mtp) 



-R-l/ 2 S + i/ 2 e 



-i(wt—rrvp) 



(2.9) 



where M+i/2 and R_i/ 2 are functions of the coordinate r, and S+i/2 and S-1/2 are functions of 
the coordinate 9. Following the steps in Refs. 0, 29], after the tedious calculation we find 



that the Dirac equations can be reduced to the following radial and angular parts 



/2) 



= (A + lj ur)VAR +1/2 
APj(VA]R +1/ 2) = (A - i/xr)R_i 

^l/2'S'+l/2 = — A5-1/2, 
Ao + i/2, 



1/2 



with 



V n ~- 
P f = 

n 

r - 

*-"n 
if = 



(2.10) 
(2.11) 
(2.12) 
(2.13) 



<9r 



zif n dA 



d iK n dA 

dr A A dr ' 



d m 

09 sin 9 
d 



m 



89 sin 9 



+ n cot 9, 



+ n cot 9, 



r 2 cu. 



We can eliminate S+1/2 (or S-1/2) from Eqs. ()2.12|) and ()2.13|) and obtain 
1 d f . „ d\ m 2 + 2mscos9 + s 2 cos 2 9 



sm9d9 \ Sine d9) 



sin 2 9 



+ s + A s 



S, = 0, 



(2.14) 



(2.15) 



where A+1/2 = A 2 — 1 and A-x/2 — A 2 . The angular equation (|2.15|) can be solved exactly and 
A s = (I — s)(l + s + 1), where I is the quantum number characterizing the angular distribution. 
Thus, for both cases s = ±1/2 we have 

A 2 =(Z + iY. (2.16) 
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In what follows, we focus our attention on the massless Dirac field. Then, we can eliminate 



i-i/2 (or 
M-l/2 (or 



-1/2) from Eqs. (j2.1(J|) and (|2.11j) to obtain a radial decoupled Dirac equation for 
-1/2)5 and we find that both them can be expressed as 



dr V dr 



+ PR S = 0, 



with 



P 



K 2 



dr 



A 



+ Asiur + — I s 



Introducing an usual tortoise coordinate 



2) ~d^ 



(2.17) 



f2.18) 



dr* = (r 2 /A)dr, 



(2.19) 



and resolving the equations in the form 



A -s/2 



(2.20) 



we obtain 



d^s 
dr 2 



+ (w 2 - V)^ s = 0, 



where 



A d 

4r 2 dr 



d /A 



dr V r 



+ 



s 2 r A 



4 



d /A 

dr \ r 4 



+ zsu;r z 



d (A 



dr \ r 



+ 



A 2 A 



(2.2i; 



f2.22) 



The potential V is complex and the complex frequency u is not separated from the potential. 
As many authors 

Q H Q 

have found that such complex potential can also present 
correct quasinormal frequencies provided we use proper boundary conditions. 



III. NUMERICAL APPROACH 

The QNMs of AdS black holes are usually defined as the solutions of the relevant wave 
equations characterized by purely ingoing waves at the black hole event horizon and vanishing 
of the perturbation at the radial infinity. Therefore, the boundary conditions on wave function 
\l/ s (or M. s ) at the event horizon (r = r+) and the radial infinity (r — > 00) can be mathematically 
expressed as 

[ A- s / 2 e~ iujr * r -> r+, 
^ s ~ A s/2 R S ~ I + (3.1) 

r — > 00. 
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Equations (|2.2ip . (|2.22|) and ()3.1|) determine an eigenvalue problem for the quasinormal fre- 
quency u of the Dirac field perturbations. 

We will calculate the quasinormal frequencies for outgoing Dirac field (i.e., for the case 
s = —1/2 [3^) by using the Horowitz-Hubeny approach 11 1. Writing $ s for a generic wave 
function as 

$ s = ^e^\ (3.2) 

then, Eq. (|2.21|) can be rewritten as 

^ (4 (3-3. 

To map the entire region of interest, r + < r < +00, into a finite parameter range, we change 
variable to x = 1/r. Define a new function B(x) as 



= X 2 _ i + + V ^ + x 3 + g2 x 4 + ^ (3 4) 



then / can be rescaled as 



2 . 



X 

In terms of the new variable x, Eq. ()3.3|) can be expressed as 



n , J 2 $ s Tlx) d$ s U(x) , , 

^tt + ^^-7^ + ( 1 \ 2 *s = 0, (3.6) 

(XX X — X-L- (XX \X — X-L- I 



where 



2 



n, \ B i X ) f 1 + X + 2 X 1 ^2 S 

S(x) = - — = ^x 2 + ^ + Q 2 x 3 

yJb — X_|_ J y X_|_ X_|_ X_|_ 

l+ad 9 £ 1 ^ 9 ,\ f3(l + xi + Q 2 xi) 9 ^ 9 
^x 2 + -=- + Q 2 x 3 - + v + x 2 -2x- 4QV - 2iu 

B(x)d 2 B(x) s 2 fdB(x)\ 2 . dB(x) x2 



£/(*) = -V = -^-^ - T ^ J + to-^ - (3.7) 

We can expand S'(x) about the event horizon x — x+ as S'(x) — X^n=o ^n(. x — x_|_) n since S(x) 
is a polynomial of degree 6, and similarly for T(x) and U(x). 

In order to evaluate the quasinormal frequencies by using the Horowitz-Hubeny method, we 
need to expand the solution to the wave function $ s around x + , 

00 

$, = {x- x + ) a a k( x - x + ) k . (3.8) 

k=0 
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It is easy to show that the index a has two solutions a = — § and a = | + ^. Because we need 
only the ingoing modes near the event horizon, that is to say, $ s must satisfy the boundary 
condition (jHHJ), we take a = — §. Then $ s becomes 



oc 



= (x - x + )- s ' 2 £ - x+) fc . (3.9) 

A:=0 

Substituting Eq. ()3.9j) into Eq. ()3.6|) we find the following recursion relation for a n 

j 71— 1 

«n = ^ [ffe - (/c - | - l) S' n _ fc + (/c - T n _ fc + U n _ k a k , (3.10) 



where Z n = 4nK[(n — s)k — iu]. 

The boundary condition ()3.1|) at r — >• oo becomes 



oo 



^a fc (-x + ) fc = 0. (3.11) 

k=0 

Now the computation of quasinormal frequencies is reduced to find roots of the Eq. (|3.11|) . 
We truncate the sum ()3.11|) at some large k = N and then check that for greater k the roots 
converge to some true roots, i.e., quasinormal frequencies. In the series f!3.10|) and (|3.11|) each 
next term depends on all the preceding terms through the recursion relations, and the roots of 
(13.11)) suffer a sharp change for a small change on any of the input parameters, especially for 
finding the higher overtone modes. In order to avoid the "noisy" we increase the precision of 
all the input data and retain 70-digital precision in all the intermediate process. 



IV. NUMERICAL RESULTS 

In this section we represent the numerical results obtained by using the numerical procedure 
just outlined in the previous section. The results will be organized into four subsections: the 
fundamental QNMs, dependence of the black hole charge, higher modes, and dependence of 
the angular quantum number. 



A. Fundamental Quasinormal Modes 

The fundamental quasinormal frequencies (n = 0) corresponding to A = 1 Dirac per- 
turbation for RNAdS black holes (r + = 5, 10, 25, 50, 75, 100) with the charge Q = 
0, 0.1Q ex t, 0.3Q ext , 0.5Q ext are given by table (JTJ) and corresponding results are drawn in 
Fig. (^). From the table and figures we find that, for the large black holes, both the real and 
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TABLE I: The fundamental quasinormal frequencies (n = 0) corresponding to A = 1 Dirac 
perturbation of the RNAdS black holes (r+ = 5, 10, 25, 50, 75, 100) with the charge Q = 
0, 0.lQ ext , O.SQext, 0.5Q ex t- 



r + 


u(Q = 0) 


w(Q = 0.1Q ext ) 


lo(Q = 0.3Q ext ) 


L0(Q = O.SQext) 


100 


199.769-152.016i 


198.415-153.628i 


187.351-169.030i 


182.694-217.611i 


75 


149.834-114.012i 


148.817-115. 220i 


140.519-126.770i 


137.024-163. 204i 


50 


99.8993-76.0043i 


99.2226-76. 8101i 


93.6914-84.5094i 


91.3576-108.796i 


25 


49.9795-37.9952i 


49.6411-38.39681 


46.8775-42.2439i 


45.7003-54.3789i 


10 


20.0744-15. 1751i 


19.9398-15.33571 


18.8398-16.86751 


18.3406-21. 6987i 


5 


10.1834-7.54760i 


10.1173-7.62708i 


9.57673-8.38092i 


9.27945-10.75511 



the imaginary parts of the quasinormal frequencies are the linear functions of the Hawking 
temperature, and the lines are described by 



Re(uo(Q = 0)) = 8.367T, 
Re(u(Q = 0.1Q ext )) = 8.309T, 
Re{cu{Q = 0.3Q ext )) = 7.845T, 
Re(uj(Q = 0.5Q ext )) = 7.652T, 



(4.1) 



and 



-Im(uj(Q = 0)) = 6.371T. 
-Im(uj(Q = 0.lQ ext )) = 6.439T, 
-Im(u(Q = 0.3Q ext )) = 7.084T, 
-Im{uj{Q = 0.5Q ext )) = 9.121T, 



(4.2) 



which show that the slope of the lines for the real parts decreases as the charge increases, but 
the slope of the lines for the magnitude of the imaginary parts increases as the charge increases. 

According to the AdS/CFT correspondence, the decay of the Dirac perturbation can be 
translated into a time scale for the approach to the thermal equilibrium in CFT. The time 
scale is simply given by the imaginary part of the lowest quasinormal frequency r = l/|7m(o;)|. 
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5 10 15 20 25 30 5 10 15 20 25 30 



FIG. 1: Graphs of fundamental quasinormal frequencies u versus the Hawking temperature T for 
black holes (r+ = 5, 10, 25, 50, 75, 100) with A = 1. The left figure is drawn for Re{oj) in which the 
lines from the top to the bottom correspond to Q = 0, 0.1Q ex t, 0.3Q ex t, 0.5Q e xti and the right one for 
Im(uj) in which the lines from the top to the bottom correspond to Q = 0.5Q ex t, 0.3Q ex t, 0.1Q ex t, 0. 
The figures show that both the real and the imaginary parts of the frequencies are the linear functions 
of T. 

From Eqs. ()4.2j) we know that for three-dimensional CFT the time scales are 

r (Q = 0) = 0.1570/T. 
t{Q = 0.1Q ext ) = 0.1553/T, 
t(Q = 0.3Q ext ) = 0.1412/T, 

r(Q = 0.5Q ext ) = 0.1096/T, (4.3) 
which show that different black hole charge presents different time scale. 

B. Dependence on the black hole charge 

The relation between the quasinormal frequencies and the charges of the black hole for 
n = 0, 1, 3, 5, r + = 100 and A = 1 is shown by Fig. (J2J). We know from figures for Re(u) with 
Q/Qext hi Fig. (EJ that the local minima of the "wiggle" appear as the scalar and gravitational 
perturbations [16] [34]. The "wiggle" becomes shallower as the overtone number n increases and 
disappears for n > 5. However, we should point out that the curves of Re(u) versus Q/Q ex t 
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; Re(w) 




Qc/Qext=0.51 5 
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FIG. 2: Graphs of to with Q/Q ex t for n = 0, 1, 3, 5, r + = 100 and A = 1. The left four figures are drawn 
for Re(u>) which show the "wiggle behavior". The right figures for Im(tv) in which the continuous 
lines describe the oscillatory modes and the dashed lines indicate the non-oscillatory modes. Values 
of Qc/Qext in the graphs indicate when non-oscillating dominates. 

'or the Dirac field are different from those of the scalar field and gravitational perturbations 

3Q. 

The figures for Im{uj) versus Q/Q ex t in Fig. Q tell us that there are two stable classes 
of solutions for QNMs when the value of the charge nears Q c , i. e., the oscillatory modes 
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; Re(w(n))/r+ 
lm(w(n))/r+ 












: Q=0.5Qext ^^^^^ 
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FIG. 3: Graphs of versus n corresponding to A = 1 Dirac perturbations of a large RNAdS black 
hole (r + = 100) for Q = 0, 0.1Q ext , 0.2Q ext , 0.3Q ex t, 0.5Q ex t, 0.6Q ex t- The dashed lines describe 
Re(u>(n)) solid lines are Im ^\ n >) _ it i s interesting to note that the quasinormal frequencies 

become evenly spaced for high overtone number n. 



(the frequencies have both the real and imaginary parts) and the non-oscillatory modes (the 
frequencies are purely imaginary). For the RNAdS black holes with the charge Q < Q c , the 
Dirac perturbation decay is dominated by oscillatory modes. But for the black holes with the 
charge Q > Q c , the decay is dominated by non-oscillatory modes. The figures also show us 
that the value of the critical value of the charge Q c slightly increases as the overtone number 
n increases. 



C. Higher modes 



In Fig. © we graph versus n corresponding to A = 1 Dirac perturbations of a large 



r+ 



RNAdS black hole (r + = 100) for Q = 0, 0.1Q ext , 0.2Q ext , 0.3Q ext} 0.5Q ext , 0.QQ ext . The 



11 



behavior of the QNMs can be written as 

u)(n) 



uj(n) 



1.35 - 2.25i)n + (2.136 - 1.50i) for Q=0, 

w(n) 



^f - ~ (1.325 -2.28i)n+ (2.12 -1.52i) for Q = 0.1Q ext , ^ ~ (1.224- 2.42i)n+ (2.14- 1.30i) 

1.8n + 1.51)z for 



for Q = 0.2Q ext , ^ ~ (1.168 - 2.56i)n + (1.90 - 1.675i) for Q = 0.3Q ext , 2g2 ^ (1.H6 



o.sg 



exti 



u)(n) 



r+ 



2.98i)n + (2.1 - 2.175i) and -(3.35n + 2.198)i for Q 
Q = 0.6Q ext . 

Cardoso et al [19j found that the scalar quasinormal frequencies for the large Schwarzschild- 
AdS black hole become evenly spaced for high overtone number n. Our results extend their 
argument to the Dirac quasinormal frequencies of the RNAdS black hole for fixed value of the 
charge. The spacings between frequencies are given by 



LO 


[n+1)- 




[n] 




r + 






LO 


[n + 1)- 




[n] 




r+ 






LO 


> + l)- 


LO 


[n) 




r+ 






LO 


> + l)- 




[n] 




r+ 






LO 




LO 


[n) 


r + 



4.35 -2.25i) 

1.325 - 2.28i) 

1.224 - 2.42i) 

1.168 — 2.56i) 

1.116 -2.98i) 
- 3.35i) 



for Q = 0, 

for Q = 0.1Q ex t, 

for Q = 0.2Q ext , 

for Q = 0.3Q ext , 

for Q = 0.5Q ext . 



(4.4) 



From which we find that, for lowly charged RNAdS black hole, the real part in the spacing ex- 
pression decreases as the charge increases, while the magnitude of the imaginary part increases 
as the charge increases. 

It is obvious that there are two stable classes of solutions for the QNMs with Q = 0.5Q ext . 
In the first class the frequencies have both the real and imaginary parts (oscillatory modes), 
and in the second class the frequencies are purely imaginary (non-oscillatory modes). But there 
is only an oscillatory mode for the case with the small charge (say Q = 0.1Q ext ), and there is 
only a non-oscillatory modes for the case with the large charge (say Q = 0.6Q ext ). 



D. Dependence on the angular quantum number 

We also study the relation between the quasinormal frequencies and the angular quantum 
number / and the results are presented in figure (pjj. We learn from the figures that the 
Re(uj) increases (decreases the oscillatory time scale) and \Im{u)\ decreases (increases the 
damping time scale) as the angular quantum number I increases for fixed n, though quasinormal 
frequencies depend very weakly on I. Besides, the figures tell us that for the imaginary parts 
\dW\ increases as the charge Q increases. We also find that the larger the overtone number n, 
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FIG. 4: Graphs of the increment dW with A (or I) for several value of n = 0, 2, 4, 6, 8 and Q = 
0-lQext, 0-3Qext- The left two figures are drawn for dW = Re(oo(l,n)) — Re(uj(0,n)) and the right 
two are for dW = —Im(uj(l,n)) + Im(oj(0, n)). The figures show that Re(u) increases and \Im(uj)\ 
decreases with increase of / for fixed n, though the dependence is very weak. Besides, the figures tell 
us that the larger the value of the charge Q, the larger the increment \dW\ for the imaginary parts. 
We also find that the larger the overtone number n, the smaller the increment \dW\ except for the 
case of Im{u){l,n = 0)) with Q = 0.3Q ex t. 
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the smaller the increment \dW\ except for the case of Im(uj(l,n = 0)) with Q = 0.3Q ext . 
V. SUMMARY 

The Dirac fields in the RNAdS black hole spacetime are separated by means of the Newman- 
Penrose formulism. Then, the quasinormal frequencies corresponding to the Dirac field per- 
turbations are studied by using the Horowitz-Hubeny approach and the results are presented 
by table and figures. We have shown that: (i) Both the real and the imaginary parts of the 
fundamental quasinormal frequency for the RNAdS black holes (r+ = 5, 10, 25, 50, 75, 100) 
with the charge Q — 0, 0.1Q ext , 0.3Q ext , 0.5Q ext are the linear functions of the Hawking tem- 
perature. The slope of the lines for the real parts decreases as the charge increases, but the 
slope for the magnitude of the imaginary parts increases as the charge increases. According 
to the AdS/CFT correspondence, the decay of the Dirac perturbation can be translated into 
a time scale for the approach to the thermal equilibrium in CFT. The time scale is simply 
given by the imaginary part of the lowest quasinormal frequency, r = 1/| Jra(u;)|. Thus, for 
three-dimensional CFT the time scales are r = 0.1570/T, 0.1553/T, 0.1412/T, 0.1096/T for 
Q = 0, 0.1Q ext , 0.3Q ext , 0.5Q ext , respectively, which show that different charge presents dif- 
ferent time scale, (ii) There are two stable classes of QNMs when the value of the charge nears 
Q c , i. e., the oscillatory modes and the non-oscillatory modes. The Dirac perturbation decay 
is dominated by the oscillatory modes for the black hole with the charge Q < Q c , but the 
decay is dominated by the non-oscillatory modes for the black hole with the charge Q > Q c . 
The value of the critical value of the charge Q c slightly increases as the overtone number n 
increases, (iii) The Dirac quasinormal frequencies of the large RNAdS black holes become 
evenly spaced for high overtone number n. For lowly charged RNAdS black hole, the real part 
in the spacing expression decreases as the value of the charge increases, while the magnitude of 
the imaginary part increases as the charge increases, (iv) The Re(u) increases (decreases the 
oscillatory time scale) and \Im(uo) \ decreases (increases the damping time scale) as the angular 
quantum number I increases for fixed n, though quasinormal frequencies depend very weakly 
on I. 
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